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We study the dynamics of a single spin embedded in the tunnehng barrier between two supercon- 
ductors. As a consequence of pair correlations in the superconducting state, the spin displays rich 
and unusual dynamics. To properly describe the time evolution of the spin we derive the effective 
Keldysh action for the spin. The superconducting correlations lead to an effective spin action, which 
is non-local in time, leading to unconventional precession. We further illustrate how the current is 
modulated by this novel spin dynamics. 

INTRODUCTION 

The analysis of spins embedded in Josephson junctions has had a long and rich history. Early on, Kulik l] argued 
that spin flip processes in tunnel barriers reduce the critical Josephson current as compared to the Ambegaokar- 
Baratoff limit Q. More than a decade later, Bulaevskii et al. P| conjectured that 7r-junctions may be formed if spin 
flip processes dominate. The competition between the Kondo effect and the superconductivity was elucidated in jj]. 
Transport properties formed the central core of these and many other pioneering works, while spin dynamics was 
relegated to a relatively trivial secondary role. In the current article, we report on new non-stationary spin dynamics 
and illustrate that the spin is affected by the Josephson current. As a consequence of the Josephson current, spins 
exhibit novel non-planar precessions while subject to the external magnetic fleld. A spin in a magnetic fleld exhibits 
circular Larmor precession about the direction of the fleld. As we report here, when the spin is further embedded 
between two superconducting leads, new out-of-plane longitudinal motion, much alike that displayed by a mechanical 
top, will arise. We term this new effect the Josephson nutation. We further outline how transport is, in turn, 
modulated by this rather unusual spin dynamics. Our predictions are within experimental reach, and we propose a 
detection scheme. 

THE SYSTEM 

The system under consideration is illustrated in Fig. ^ It consists of two identical ideal s-wave superconducting 
leads coupled each to a single spin; the entire system is further subject to a weak external magnetic fleld. In Fig.l^, 
A*L,fl denote the chemical potentials of the left and right leads, B is a weak external magnetic fleld along the z-axis, 
and S = {Sx, Sy, Sz) is the operator of the localized spin. The Hamiltonian of the system reads 

H = Ho+Ht, Ti-o — Hl +^-11 - ^BzSz, (1) 

Ht^ ^fl.ka [ToSaa, + Ti Cr„„, • S ] CLpa' + h.C, (2) 

k,p,a,a' 

where Ti^ and Tin are the Hamiltonians in the left and right superconducting leads, while c]^^^ (cika) creates (an- 
nihilates) an electron in the lead a in the state k with spin a in the right/left lead for i = L/R respectively. The 
. ^ ' vector cr represents the three Pauli matrices and /i is the magnetic moment of the spin. When a spin is embedded in 
, the tunneling barrier, the conduction electron tunneling matrix becomes, not too surprisingly, spin-dependent 0, Q 
5^ • T — [TqI + TiS-o-]. Here Tq is a spin- independent tunneling matrix element and Ti is a spin-dependent matrix element 



FIG. 1: Magnetic spin coupled to two superconducting leads. 
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originating from the direct exchange coupling J of the conduction electron to the localized spin S. We take both 
tunneling matrix elements to be momentum independent. This is not a crucial assumption and is merely introduced 
to simplify notations. Typically, from the expansion of the work function for tunneling, |f ~ J/U ^ where U is the 
height of a spin-independent tunneling barrier A weak external magnetic field ^ IQQGauss will not influence 
the superconductors and we may ignore its effect on the leads. In what follows, we abbreviate by B. The operator 
gZ0/2 jg ^j^g (single electron) number operator. When the junction is linked to an external environment, the coupling 
between the junction and the environment induces fluctuation of the superconducting phase (0). 



THE EFFECTIVE ACTION 



Josephson junctions are necessarily embedded into external electrical circuits. This implies that the dynamics will 
explicitly depend on the superconducting phase (j). The evolution operator is given by the real-time path integral 



Z = j D(j)DS exp [i5circuit(0) + iSsphiiS) + i5tunnol(0, S)] 



(3) 



The effective action 5tunnci describes the junction itself. We generalize the formerly known effective tunneling action 
for a spin-less junction ^ J^j to the spin-dependent arena to obtain 



^tunnel = -2 i dt j, dt'a{t,t') [T^ + TfS{t) ■ S{t')] cos 
Jk Jk 

- 2(f dt(f dt' f3{t,t') [T^ - T^S{t) ■ S{t')] cos 
Jk Jk 



m + m 



IK JK 

where ia{t,t') = G{t,t')G{t' ,t), i(3{t,t') = F{t,t')F\t,t') and the Green functions 

G{t,t')^-^Y,{TKC^^{t)cUt')) 

k 

F{t,t') = -iY,{TKC^.^{t)c^^.iit')) 

k 

FHt,t')^-zY,{TKci,{t)cl^,{t')) . 



(4) 

(5) 
(6) 
(7) 



In Eq. §^ denotes integration along the Keldysh contour. We now express the spin action on Keldysh contour 
in the basis of coherent states 
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spin — f diB • S + SwZNW 
K 



(8) 



Here, S denotes the magnitude of the spin S. The second, Wess-Zumino-Novikov-Witten (WZNW), term in E 
depicts the Berry phase accumulated by the spin as a result of motion of the spin on a sphere of radius 
Explicitly, 

SwzNW j( dt[S{t,T) ■ (9rS(t,r) x dtS{t,T))] . 



(9) 



The additional integral over r allows us to express the action in a local form. At r = the spin is set along the z 
direction at all times, S{t, 0) = const; at r = 1 the spin field corresponds to the physical configurations, S(t, 1) = S{t). 



DYNAMICS 



We now perform the Keldysh rotation, defining the values of the spin and the phase variables on the for- 
ward/backward branches of the Keldysh contour (e.g., S"'' for the upper and lower branch) and rewriting all the 
expressions in terms of their average (classical component S) and difference (quantum component 1): 



S = (S'' + S')/2, l^S^'-S', S-1 = 0. 



(10) 
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FIG. 2: The sphere or radius S for the vectors S"'\t) is shown. The path C describes the evolution of the spin along the upper 
(m) and lower (l) branches of the Keldysh contour. To properly describe the spin dynamics on this closed contour we analyze 
the WZNW action, see Eq. @. For clarity, we draw a small piece of the closed trajectories. 



After the Keldysh rotation we obtain 

SwzNW^^J^ dTy"di[S"(i,r)-(a,S"(t,r) xatS"(i,r))-(«^0] • (h) 

The relative minus sign stems from the backward time ordering on the return part of C. The individual WZNW 
phases for the upper (u) and lower (l) branches are given by the areas spanned by the trajectories S"''(i) on the 
sphere or radius S divided by the spin magnitude (S). The WZNW term contains odd powers of 1. Insofar as the 
WZNW term of Ea. pi|l is concerned, the standard Keldysh transformation to the two classical and quantum fields, 
S and 1, mirrors the decomposition of the spin in an anti-ferromagnet (AF) to the two orthogonal slow and fast fields 
[l5| . The difference between the two individual WZNW terms in Ea. (|ll|l is the area spanned between the forward 
and backward trajectories. For close forward and backward trajectories the WZNW action on the Keldysh loop may 
be expressed as 



^ J dt\-iSxdtS). (12) 



SwZNW — g2 

For the spin part of the (semi-classical) action we, then, obtain 

5.pi„ = J dtB-l+-^ J dtl-{SxdtS). (13) 

Next, we perform the Keldysh rotation to the classical and quantum components with respect to both the phase and 
spin variables in the tunneling part of the effective action. Towards this end, we introduce (with notations following 
Refs. |l[l3) 

0=(0"+/)/2 , x = • (14) 

With these definitions in hand, the tunnehng part of the action reads 

'^tunnel — 

Sc+Sp , (15) 

where the normal (quasi-particle) tunneling part Sa is expressed via the Green functions = 9{t — t'){a^ — a^) and 
a^{uj) = a> +a<, where ia>{t,t') = G> {t,t')G< {f ,t) and ia<{t,t') = G<{t,t')G> {t\t). Similarly the Josephson- 
tunnehng part 5/3 is expressed via the off-diagonal Green's functions = 9{t — t'){j3^ — /?<) and fi^ [ijj) = /3> + /3<, 
where il3>(t,t') = F> {t,t')F^> {t,t') and il3<{t,t') = F<{t,t')F^<(t,t'). In this paper we are interested in the 
interaction between the supercurrent and the spin. Thus we provide the expression for the Josephson part: 



dt J dt' Afi^{t,t') X 

{ [2Ti - 2T?S(,) . S(01 ..n ^ COS - i T? 1(0 ■ 1(0 cos ^21 xgOj M±*« 
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jrf m ■ S{t') COS ^ COS ^ - Ti Sit) ■ m sin ^ sin ^| cos ■ 
+ I dt [ dt' P^it,t') X 



I [4T^ - 4rf S(t) . Sit')] sin ^ sin ^ + • l(t') cos ^ cos ^| 



+ ^it') 

cos - 



2rf lit) ■ Sit') COS ^ sin 41 + 2Tf S(0 • l(t') sin ^ cos ^ ^^'^ + ^^^'^ 



2 

(16) 



The normal-tunneling part Sa is obtained from S/b by the following substitution: /3^/^(<, i') a^^^ it,t'), 4>it') 
—(pit'), and x(^') ~^ ~xit')- The Keldysh terms (those including and a''^), which normally give rise to random 
Langevin terms (see, e.g., Ref. |lOj ) are, in our case, suppressed at temperatures much lower than the superconducting 
gap (r <C A), due to the exponential suppression of the correlators P'^ iw) and iui) at < A. 
To obtain we start from the Gor'kov Green functions 

F>(i,t') = -z5:^e-^^-(*~*') , F>t(^,^')-^E^^■^''^'*"*'^ (l^) 

k k 



where the quasi-particle energy Ek = \/ A^ + e^, being the free-conduction-electron dispersion in the leads. Putting 
all of the pieces together, we find that 

a^it - t') = eit - t') J2 sin HE, + E,)it ~ t')] . (18) 

The kernel /3^(t — t') decays on (short) time scales of order Oih/A). Varying the full action with respect to the 
quantum components 1 and x ^ind setting these to zero, we obtain coupled equations of motion for both the spin and 
phase: 



= Sit) xB + T^ J dt' 4/3^(i - t') Sit) X Sit') cos '^^^^ ^ '^^^'^ , (19) 

1x^0 = -Jdt' 2p^it - t') iTi - T^Sit) ■ Sit')) sin . (20) 

Note, that, if the rest of the circuit contains dissipative elements, e.g., resistors, then ^circuit will contain the non- 
vanishing Keldysh components, and one should include the corresponding Langevin terms into Eq. (|2()|l . The rather 
complicated equations of motion (|19|l and (|2()|l are very general. To make headway, we now adopt a perturbative 
strategy. In Ea. (|19(l . we first assume an ideal voltage bias, i.e., an imposed phase 0(i) = ujjt, where the "Josephson 
frequency" uij = 2eV/fi. To this lowest order, we neglect the influence of the spin on the phase. Next, we use the 
separation of characteristic time scales to our advantage. To this end, we note that the spin dynamics is much slower as 
compared to electronic processes, i.e. wj,i?<C A. This separation of scales allows us to set S(i') ~ Sit) + it' — t)dS/dt 
in the integrand of Eq. (jl9|l , wherein we obtain 

dS dS 

— = AS X — sinwjt-l-S X B. (21) 

(XC civ 



Here, 



00 cc 

-4r2 J dt-t- P^'it) • sin w -2T?L0j j dt-t'- p"it) = 2wj ^ 



lAplTil 



EkEpiEk+Ep)3 

u k,p ^ 

- 9T2^2^ f f dzidz2 _ TT^ ^2 „2 _ ffl fr,r,^ 

" 1^ A J J (coshzi-f coshz2)3 " 8 1^ A " 32 A ' ^ ' 

with gi = (27rTip)^ the spin channel conductance. In Ea. (|22|l we employ the separation of time scales (tjj <C A) 
again. When expressed in the spherical coordinates (in the semi-classical limit) S — S'(sin0 cos 0, sinfl sin cos 0), 
Ea. (|21|l transforms into two simple first order differential equations 

^ = ^ (23) 

dt 1 + S^X^sin^iujjt) ' ^ ' 

dO ^d(t> , ^ , 

— = -5'A-y-sin6lsinu;ji . (24) 

(Jill lLv 



FIG. 3: The resulting spin motion on the unit sphere in the general case. As in the motion of classical spinning top, the spin 
exhibits undulations along the polar direction. 



These equations can be solved exactly. For a spin oriented at time t = at an angle 9o relative to B, 

B 



m = - 



: tan- 1 [^/Y+S^i{m{uJJt)] , 



e{t) = 2 tan 



(1 - ccos(wjt))(l + cy 



(1 + ccos(wjt))(l — c) 



tan — 
2 



withc= S'A/VT+S^ and 7 = -SXB/2ujjc. For 5A < 1, ~ -Bt and 9- On-- SX{B/ uj) sin 0n cos ojjt. Typically, 
whenever a spin is subjected to a uniform magnetic field, the spin azimuthally precesses with the Larmor frequency 
ujl — B. In a Josephson junction, however, the spin exhibits additional polar (9) displacements. The resulting 
dynamics may be likened to that of a rotating rigid top. The Josephson current leads to a non-planar gyroscopic 
motion {Josephson nutations) of the spin much like that generated by torques applied to a mechanical top. For small 
A, we find nutations (see Fig.O of amplitude 6i — 02 (x sin 6 ■ S ■ X ■ (x sin6 ■ S ■ gi ■ ^. 

The origin of the first term on the rhs of Eq. (|21|l can be understood as follows (this origin can be also traced in the 
calculations): the spin is subject to the electron-induced fluctuating field h = TiY.e^'^^^c' crc + h.c. The same coupling 
may be thought of as an influence of the spin on the leads, which results in a non-zero low-frequency contribution 
(5hto h. Since the response function of the electron liquid is isotropic but retarded, 5h(t) is not aligned with S{t) but 
contains information about the values of S{t') at earlier times. The response function decays on a time scale ~ ?i/A, 
much shorter than the period of the spin precession, ^ 1/B. As a consequence, in addition to a contribution cx S the 
field h acquires a component oc S / A, which leads to the first term on the rhs of Eq. (|21|l . 

The rhs of the second equation of motion H2U|) clearly corresponds to the Josephson current. Indeed, in the Keldysh 
formalism one has / = (27r/$o)95/9x (instead of / = {2Tr/^o)dS /d(f>). Thus we obtain for the Josephson current 



27r 

$0 



dt' 213" {t - t') {T^ ~ TfS{t) ■ S{t')) 



m + m 



We start from the lowest order (local in time) adiabatic approximation, i.e. we set S(i) = S(i') and (/)(t) = 0(t') 
yields 



. / N 27r 
5>0 



rp2 



sin ip{t) 



(25) 



This 



(26) 



where Ej^ = 2Tg J dtf3"{t) = n'^p'^T^A = (l/4)goA is the spin-independent Josephson energy [2| {go being the 
conductance of the spin- independent channel). The second term of Eq. H26I) gives the spin-related reduction of the 
Josephson critical current studied in Ref. We now evaluate the lowest-order correction to this equation due to 
deviations from locality in time and spin precessions. Expanding S{t') in Eq. H25I) in {f — t) and using the fact that 
for the Larmor precession we have S • S = and S ■ S = B^ {Si — S^), we find a correction to the Josephson current 
which depends on S^: 



. / ^ 27r 
5>o 



E 



J,0 



rp2 



6Ej{Sl-S^) 



sin(/){t) , 



(27) 



6 




FIG. 4: A SQUID-based detection scheme. The SQUID monitors the magnetic field produced by the magnetic cluster in one 
of the junctions. 

where SEj = -T^B^ J dtf3^{t)t^ = {7r^/16)T^p^ (^VA). Here, we clearly elucidated the manner in which the spin 
dynamics alters the Josephson current. 

For S =1/2 the semi-classical approximation is insufficient. In this case it is easier to perform a calculation with spin 
operators Il3!|, rather than a path integral. One, then, obtains |l3l| an expression for the Josephson current identical 
to Eq. H25() with S{t) being, however, the spin operator in the interaction representation. Using the commutation 
relations of the spin operators one obtains an extra contribution to the Josephson current proportional to Sz- This 
allows reading out of the spin state via the Josephson current. This extra contribution scales as S while the spin 
dependent contributions in Eq. l|77| scale as S"^. 

DETECTION 

We now briefly discuss a detection scheme for the Josephson nutations for S ^ 1, e.g., in the semi-classical limit. In 
principle the nutations should affect the Josephson current. The level of approximation employed in this paper was, 
however, insufficient to describe this effect. Indeed, one has to substitute S(t) containing the nutations into Eq. H25|l . 
As the amplitude of the nutations is of the order gi, the correction to the current will be of the order gf. We will 
study this correction elsewhere. 

Here we discuss a more direct detection strategy. The spin motion generates a time-dependent magnetic field, 
(5B(r,t) = |^[3r(r-m(t)) — r^m(t)]/r^, superimposed on the constant external field B. Here, r is the position relative 
to the spin, with magnetic moment m{t) — /iS(t). A ferromagnetic cluster of spin S = 100 generates a detectable 
field SB ~ 10~^° T appears a micron away from the spin. For a SQUID loop of micron dimensions located at 
that position, the corresponding flux variation 5$ ^ 10~^<I>o (with $o a flux quantum) are within reach of modern 
squid's. For such a setup with Ti/Tq ^ 0.1, the typical critical Josephson current is Jg^^ ~ 10 pA, |A| = 1 meV, 
and eV ~ 10~'^|A|, we find that AS* ^ 0.1. Since Sx = S'sin^cos^, Sy = S sinO sm(f), the spin components orthogonal 
to B vary, to ffist order in (AS'), with Fourier components at frequencies \lul ii^jj {^^l = B), leading to a discernible 
signal in the magnetic field B + (5B. For a field B ^ 200 G, ujl ^ 560 MHz, and a new side band will appear at 
\ujl — wjI, whose magnitude may be tuned to 10-100 MHz. This measurable frequency is markedly different from the 
Larmor frequency col- 

The efficiency of the detector may be further improved by embedding the spin in one of the Josephson junctions 
of the SQUID itself. The setup is sketched in Fig. 0] The Josephson junction containing the spin is used both for 
driving the nutations and, together with the second junction of the SQUID, for detecting them. 
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DISCUSSION 



In this article, we illustrated that the dynamics of a spin embedded in a Josephson junction is richer than appreciated 
hitherto. We reported unusual non-planar spin motion (in a static field), which might be probed directly and which 
was further shown to influence the current in the Josephson junction. Using a path-integral formalism, we described 
this non-planar spin dynamics and the ensuing current variations that it triggers. To describe the time evolution 
we derived the effective action for a spin of arbitrary amplitude S on the Keldysh contour. In passing, we noted a 
similarity between the resultant effective action and that encountered in quantum antiferromagnetic spin chains. Our 
central results are encapsulated in the effective action Hlfifl . 

In the semi-classical limit of large S, relevant to ferromagnetic spin clusters we obtained two coupled equations 
of motion (Eqs. H19(l and H20|ll. These equations may be solved perturbatively, as outlined above, or numerically. We 
presented an exact limiting-case solution and illustrated how the new spin dynamics may be experimentally probed. 

The formalism developed can also be applied to the minimal S = 1/2 system. In this case, however, it is simpler 
to perform a calculation with spin operators , rather than a path integral. 
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